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We describe a recently discovered stable planar surface of silicon, Si(114). This high-index surface, 
oriented 19.5° away from (001) toward (111), undergoes a 2 x 1 reconstruction. We propose a 
complete model for the reconstructed surface based on scanning tunneling microscopy images and 
first-principles total-energy calculations. The structure and stability of Si(114)-(2 x 1) arises from 
a balance between surface dangling bond reduction and surface stress relief, and provides a key to 
understanding the morphology of a family of surfaces oriented between (001) and (114). 



(A preprint with high-resolution figures is at ittp: //cst-www .nrl .navy .mil/papers/sill4 .ps 



The (001) and (111) surfaces of silicon are its primary 
growth and cleavage planes, respectively. Surfaces with 
orientations between these two planes are both intrinsi- 
cally interesting and technologically important. For these 
intermediate orientations, surface morphology depends 
strongly on the orientation angle. For example, small 
misorientations away from the low-index planes produce 
vicinal surfaces consisting of terraces separated by steps 
, whereas larger angles often lead to sawtooth- like, 
grooved surfaces [f|||]. Between (001) and (111), only two 
orientations are known to form planar surfaces with sta- 
ble reconstructions: Si(113) §j| and Si(5 5 12) §. Re- 
cently, (114) surfaces have been observed both on cylin- 
drical Si samples |1 and within etch pits formed on 
Si(001) [fTpf , In addition, sawtooth-like structures com- 
posed of (114) and (113) facets have been observed on 
samples oriented between (114) and (113) ||J|. These 
results suggest that Si(114) is a planar surface that is 
thermodynamically stable against faceting. 

The (114) plane can be considered a highly vicinal 
(001) surface, oriented 0=19. 5° from (001) toward (111). 
The structure of the ideal bulk-terminated (114) surface 
is an alternating sequence of one- and two-unit-cell wide 
(001) terraces separated by double-layer steps, as shown 
in Fig. 1(a). This basic terrace-plus-step morphology also 
describes real surfaces over a large range of orientation 
angles: from slightly vicinal (001) surfaces with large 
variable- width terraces to steeper orientations consisting 
of three-unit-cell wide terraces separated by double-layer 
steps (0=11.4°) jll] . It is interesting to ask at what orien- 
tation angle this morphology breaks down, and for what 
reasons. For example, the Si(113) surface (0=25.2°) has 
a stable 3x2 reconstruction that cannot be simply de- 
scribed as the expected sequence of one-unit-cell terraces 
and double-layer steps 0. For a surface with (114) orien- 
tation, one is therefore led to consider three possibilities: 
(i) a periodic surface with a terrace-plus-step morphol- 
ogy; (ii) a periodic surface without a terrace-plus-step 
morphology, similar perhaps to Si(113)-(3 x 2); or (iii) a 
non-periodic (e.g. grooved) surface. In this Letter we use 
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FIG. 1. (a) Side view of the Si crystal lattice, showing one 
unit cell of the bulk-terminated (114) surface. Dashed lines 
emphasize the alternating single- and double- width (001) ter- 
races separated by double-layer (DL) steps, (b) Filled-state 
STM gray-scale image (210 x 180 A) of the clean Si(114) sur- 
face, (c) Atomic-resolution, filled-state image (65 x 54 A) 
highlighting the three types of observed row structures. Unit 
cells are drawn to indicate two types of local symmetry: 2x1 
[left] and c(2 x 2) [right]. The 2 x 1 unit cell is 16.3 x 7.7 A. 
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scanning tunneling microscopy (STM) to definitively es- 
tablish that Si(114)-(2 x 1) is a planar surface with a sta- 
ble terrace-plus-step reconstruction. We propose a com- 
plete structural model for the reconstructed surface, and 
use first-principles electronic-structure methods to pro- 
vide strong theoretical evidence in support of this model. 
The structure of Si(114) reveals it to be a natural exten- 
sion of vicinal (001); indeed, it marks the endpoint of a 
family of surfaces having closely related structure. 

The clean Si(114) sur face H appears in constant- 
current STM images as a well-ordered periodic array 
of row structures oriented along the [TlO] direction 
[Fig. 1(b)]. The period of this structure in the [221] di- 
rection is approximately 16 A, equal to the length of the 
bulk-terminated (114) unit cell. Along each row, the pe- 
riod in the [TlO] direction is 7.7 A, twice that of the bulk- 
terminated unit cell. Atomically resolved images reveal 
three distinct types of row structures, labeled I, II, and 
III in Fig. 1(c). The phase relationship along [TlO] be- 
tween these row structures varies over different regions 
of the surface, giving rise to areas of local 2 x 1 as well 
as c(2 x 2) periodicity [see Fig. 1(c)]. 

Our proposed structural model for Si(114)-(2 x 1) is 
closely related to the structure of vicinal Si(001). The 
Si(001) surface undergoes a reconstruction in which the 
surface atoms pair up and bond to reduce the number of 
surface dangling bonds, forming parallel rows of dimers. 
For vicinal surfaces with 9 greater than 4-5°, the great 
majority of resulting (001) terraces consist of dimers ori- 
ented parallel to the step edges (B-type terraces), with 
these B-type terraces separated by double-layer steps 
(Db steps) Q. Chadi first proposed that Db steps are 
rebonded, i.e. they incorporate another row of atoms at 
the step edge to reduce the number of dangling bonds 
|^3| . This basic morphology — B-type terraces separated 
by rebonded Db steps — has also been observed in STM 
studies of curved Si samples having regions with local ori- 
entation 9 ~ 9-11°, where the (001) terrace width is only 
three or four unit cells jll| . We propose that a variant of 
this basic morphology also describes the very narrow ter- 
races found on Si(114). Specifically, our structural model 
for Si(114)-(2 x 1), shown in Fig. 2, is an alternating se- 
quence of single- and double- width B-type (001) terraces, 
separated by rebonded Db steps (stepping up from the 
double- width terrace) and non-rebonded Db steps (step- 
ping down from the double- width terrace) . The presence 
of both rebonded and non-rebonded steps provides an im- 
portant mechanism for stress relief within the unit cell, 
as discussed in detail below. 

This structural model is strongly supported by the re- 
sults of first-principles calculations of equilibrium surface 
geometry and surface energy, and the resulting simulated 
STM images. The calculations were performed in a su- 
percell slab geometry (inversion symmetric) with ten lay- 
ers of Si and a vacuum layer equivalent to six Si layers. 
Full structural relaxation was performed using the Corn- 
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FIG. 2. Top and side views of proposed model for 
Si(114)-2(x 1) reconstruction, with fully relaxed geometry. 
Labeled structures correspond to the STM features seen in 
Figs. 1 and 3. A dimer plus two non-rebonded step-edge 
atoms is called a "tetramer." 



ing electronic-structure code of Allan, Teter, and Payne 
|l4| , which solves the Kohn-Sham equations in the local- 
density approximation (LDA) with a plane-wave basis 
and norm-conserving pseudopotentials [ I5| . The kinetic- 
energy cutoff was 10 Ry, and a single fc-point was used 
throughout. Structural relaxation was performed on all 
the atoms until the rms force was less than 0.1 eV/A Jl6| . 

The fully relaxed surface geometry for our proposed 
model includes three prominent row structures per unit 
cell, corresponding to the three types of rows seen in the 
STM images. Row I in the images arises from the line 
of rebonding atoms at the rebonded Db step edge; row 
II arises from the row of tetramers (i.e. the combina- 
tion of the dimers on the double-width terrace and the 
non-rebonded Db step-edge atoms); and row III arises 
from the dimers on the single-width terrace. The corre- 
spondence between theory and experiment is examined 
in detail in Fig. 3. Figs. 3(a) and (b) show STM im- 
ages for filled and empty states, respectively, from the 
same surface area; Figs. 3(c) and (d) show numerically 
simulated STM images (constant current surfaces) cal- 
culated by integrating the local-state density over filled 
and empty states, respectively, near the Fermi level. In 
the filled-state STM image [Fig. 3(a)], rows I, II, and III 
all show 2ao periodicity along [110] (arj=3.84 A is the 
surface lattice constant): rows I and III each appear as 
an array of maxima with 2ao spacing, whereas row II 
has a more variable appearance. In the empty-state im- 
age [Fig. 3(b)], all three rows also exhibit 2ao periodicity, 
but the features are more diffuse within each row. 
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FIG. 3. (Left) STM images (33 x 31 A) of Si(114) mea- 
sured for: (a) filled states, 1.2 V; (b) empty states, —1.7 V 
V. (Right) Simulated STM images calculated for: (c) filled 
states; (d) empty states. Black dots indicate projected atom 
positions. The simulated images are for a surface with c(2 x 2) 
symmetry. 

The general appearance of the rows is accurately re- 
produced in the simulated STM images based on our 
model. In the simulated filled-state image [Fig. 3(c)], 
row I is most prominent, and the expected 2ao periodic- 
ity is apparent in all three rows, in good agreement with 
the STM data. Along row I, weak maxima are visible be- 
tween the dominant peaks spaced 2ao apart; these max- 
ima are only occasionally resolved in STM images. Row 
II has a zigzag-like internal structure, consistent with the 
variable appearance observed experimentally. In the sim- 
ulated empty-state image [Fig. 3(d)], the 2ao periodicity 
is still evident in row II, but is much weaker in rows I 
and III, as observed in the experimental image. Finally, 
the experimental observation of areas with local 2x1 
and c(2 x 2) symmetry is also consistent with our model. 
We have calculated total energies of both reconstructions 
and find differences only at the level of 1-2 meV/A 2 . 

Several of the topographic features described above 
are a consequence of a structural distortion that occurs 
within each row. For example, at equilibrium the dimers 
that constitute row III are each tilted (buckled) so that 
the two dimer atoms have a theoretical height difference 
of Az=0.17 A. This buckling occurs on flat Si(OOl) as 
well, and is energetically favored by the transfer of an 
electron from the lower to upper atom of the dimer. A 
similar buckling occurs along the line of rebonding atoms 
that constitute row I (Az=0.18 A). Both of these buck- 
lings, although small, give rise to a large corrugation of 
the constant-current surfaces and are critical to explain- 
ing the appearance of rows I and III as lines of maxima 
with 2ao spacing. The tetramers, which constitute row II, 



have a more complicated geometry that gives rise to the 
internal structure visible in both the experimental and 
simulated images. In the equilibrium surface geometry 
shown in Fig. 2, the tetramers are essentially symmetric 
(no buckling) . In the course of relaxing from slightly dif- 
ferent initial geometries, however, we found a variation 
on this structure (also at equilibrium and with essentially 
the same total energy) with a highly buckled tetramer 
(Az=0.30 and 0.15 A for the dimer pair and step-edge 
atom pair, respectively). The simulated images in Fig. 3 
are from a surface with buckled tetramers; the similar to- 
tal energies suggest the possibility of different metastable 
geometries within row II, explaining its variable appear- 
ance in the experimental images. 

Experimentally, we observe that Si wafers oriented 
within 0.5° of (114) form large terraces of (114) sepa- 
rated by steps. Wafers with orientations a few degrees 
off (114) toward (111) form grooved surfaces consisting 
of (114) and (113) facets fl7l| , consistent with earlier x- 
ray scattering results |5| and ex- situ atomic force mi- 
croscopy (|]. These observations suggest that the sur- 
face energy of (114) is quite low, perhaps comparable 
to Si(113)-(3 x 2). We have investigated the stability 
of our Si(114)-(2 x 1) reconstruction model, relative to 
low-index reconstructed surfaces, by calculating surface 
energies, E s = [E t (N) - iV£ t bulk ]/2. Here E t (N) is the 
total energy of the supercell containing N atoms, £f ulk is 
the total energy per atom of bulk Si (computed using the 
same lattice vectors, energy cutoff, and zone sampling as 
for the surface calculations), and the factor of two ac- 
counts for the two surfaces per unit cell. We find the 
surface energies for the fully relaxed 2x1 reconstruc- 
tions of Si(lll), (001), and (114) to be 84, 86, and 85 
meV/A 2 , respectively, indicating that Si(114)-(2 x 1) is 
approximately as stable as the low-index (lll)-(2 x 1) 
and (001)-(2 x 1) surfaces. 

While our proposed model — a terraced surface with 
Db steps alternately rebonded and non-rebonded — was 
guided by our STM results, we have also carried out fully 
relaxed total-energy calculations for a number of alter- 
native structures. All lead to simulated STM images in 
strong disagreement with experiment, and most lead to 
significantly higher surface energies than our proposed 
model. Two of these models, however, shed light on 
the mechanism that ultimately stabilizes Si(114): (Model 
A) a terraced surface with every Db step rebonded; 
and (Model B) a terraced surface with every Db step 
non-rebonded. Because rebonded steps have fewer sur- 
face dangling bonds than non-rebonded steps, one might 
naively predict Model A to be strongly favored relative 
to both Model B and our model. When fully relaxed, 
however, we find that Model A has a surface energy of 
87 mc v/A 2 , very slightly higher than our model (Model 
B is significantly higher, at 97 meV/A 2 ). 

We believe that the alternation of rebonded and non- 
rebonded step edges is ultimately favored because it leads 
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to an optimal balance between surface dangling bond re- 
duction and surface stress relief. For isolated Db steps 
on Si(OOl), rebonded steps are energetically preferred 
to non-rebonded steps because the number of dangling 
bonds is reduced. However, rebonded Db steps con- 
tribute a tensile stress in the direction perpendicular to 
the step edge. For sufficiently wide terraces, this stress 
is elastically relieved, so that rebonded steps are favored. 
For very short terraces, however, this relief mechanism is 
not available, and the resulting energy penalty may be 
larger than the energy gain from dangling-bond reduc- 
tion. Si(114) balances these competing effects by incor- 
porating both a rebonded and a non-rebonded Db step in 
each unit cell, allowing for both dangling-bond reduction 
and stress relief within the unit cell. 

This mechanism of stress relief can be demonstrated 
quantitatively by calculating the LDA surface stress ten- 
sor for each model, cry = A~ 1 dE s /deij, where e$j is the 
surface strain tensor and A is the surface cell area [ ^8[ . 
Based on the above arguments, for Model A one expects 
a large tensile stress in the x-direction (i.e. perpendicular 
to the step edge) , whereas for Model B and our model one 
expects a smaller (compressive or tensile) stress. Indeed, 
we find the values of o~ xx to be 0.27 eV/A 2 (tensile) for 
Model A, -0.03 eV/A 2 (compressive) for Model B, and 
0.11 eV/A 2 for our model. These numerical results are 
consistent with the qualitative scenario described above: 
Model A has a low density of dangling bonds but high 
surface stress, which incurs a large energy penalty; Model 
B has a low surface stress but high density of dangling 
bonds, which also incurs a large energy penalty. Our pro- 
posed model compromises on both dangling bonds and 
surface stress, and by so doing is the most favorable. 

The competition between the energy gain from 
dangling-bond reduction and the energy penalty from 
tensile stress at rebonded Db steps provides a frame- 
work for understanding the surface morphology of a sur- 
prisingly large class of orientations, from slightly vicinal 
Si(001) up to Si(114) (0=19.5°). For orientations up to 
(114), the stress at rebonded Db steps can be relieved ei- 
ther on the terraces or by the presence of non-rebonded 
steps, and consequently the terrace-plus-step morphol- 
ogy is favored. For example, curved surfaces with lo- 
cal (119) and (117) orientations (0=8.9° and 11.4°, re- 
spectively) form well-ordered B-type terraces separated 
exclusively by rebonded Db step edges |IT| , |l9| . More- 
over, planar surfaces with (115) orientation (0=15.8°) 
consist of a complicated arrangement of B-type terraces 
separated by both rebonded and non-rebonded steps [JL7| . 
Beyond (114), however, the terraces become too narrow 
to afford sufficient stress relief, and so the terrace-plus- 
step morphology is no longer favored, leading to funda- 
mentally different surface morphologies. For example, 
surfaces with orientations between (114) and (113) form 
mesoscale sawtooth- like grooves [pLpl, and surfaces with 
(113) orientation are stabilized by a more complex 3x2 



reconstruction (|J^]. We therefore conclude that Si(114) 
marks the endpoint of a family of orientations, each mem- 
ber of which has a related surface morphology consisting 
of (001) terraces separated by Db steps. 
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